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We calculate the spatial entanglement between two electrons trapped in a nanostructure for a broad class 
of confinement potentials, including single and double quantum dots, and core-shell quantum dot structures. 
By using a parametrized confinement potential, we are able to switch from one structure to the others with 
continuity and to analyze how the entanglement is influenced by the changes in the confinement geometry. 
We calculate the many-body wave function by 'exact' diagonalization of the time independent Schrodinger 
equation. We discuss the relationship between the entanglement and specific cuts of the wave function, and 
show that the wave function at a single highly symmetric point could be a good indicator for the entanglement 
content of the system. We analyze the counterintuitive relationship between spatial entanglement and Coulomb 
interaction, which connects maxima (minima) of the first to minima (maxima) of the latter. We introduce a 
potential quantum phase transition which relates quantum states characterized by different spatial topology. 
Finally we show that by varying shape, range and strength of the confinement potential, it is possible to induce 
strong and rapid variations of the entanglement between the two electrons. This property may be used to tailor 
nanostructures according to the level of entanglement required by a specific application. 
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I. INTRODUCTION 

Entanglement is a quantum property which provides the 
possibiUty for quantum information/computation to overcome 
some of the limitations of traditional devices. For this rea- 
son entanglement is now considered a physical resource. Re- 
cently, semiconductor quantum dots (QDs) have been pro- 
posed as promising hardware to perform quantum informa- 
tion/computation within solid state.-^^ Their potential advan- 
tages include the existence of an industrial base for semicon- 
ductor processing and flexibility in driving the computational 
degrees of freedom by applied electro-magnetic fields and 
purposely designed trains of laser pulses.^ The system param- 
eters may be tuned, making it possible to tailor the properties 
of semiconductor nanostructures.^** The rapid technological 
advances seem to promise sophisticated engineering of QD- 
based structures, with the potential for the production of scal- 
able and coupled QD systems. A crucial requirement is then 
the possibility of generating and manipulating entanglement 
within these structures. 

In QD systems entanglement could be controlled by exter- 
nally applied electro-magnetic fields,^ or by varying nanos- 
tructure parameters. Recently the effect of the interdot 
distance on the entanglement of two electrons trapped in 
(In,Ga)As/GaAs QD molecules has been studied,^ as well as 
the effect of ionization on the entanglement of two electrons 
in a single QD.' 

In the present paper, we investigate how the geometrical 
changes in the confinement potential of single, core-shell and 
double QD structures influence the spatial entanglement"' be- 
tween two electrons trapped within the nanostructure. To this 
aim we will use the two-center power-exponential potentiaU^ 
which allows to change the confinement potential with conti- 
nuity from one structure to the others. 

We will show that small variations in the confinement po- 
tential can induce large changes in the entanglement and 
present a potential quantum phase transition between states 
with minimum and maximum entanglement. We will per- 



form a detailed study of the role of the Coulomb interaction 
in determining the entanglement and discuss the features of 
the many-body wave function which characterize the various 
entanglement regimes. 

The paper is organized as follows: the theoretical model 
is presented in section II; in section III we present and ana- 
lyze the results for the system ground state energy and related 
Coulomb interaction; section IV includes the discussion of the 
entanglement; in section V we analyze the characteristics of 
the many-body wave function; in section VI we present entan- 
glement indicators; and in section VII we discuss a potential 
quantum phase transition. Finally section VII is devoted to 
conclusions and summary. 



II. THEORETICAL MODEL 

We consider two interacting electrons confined within a 
nanostructure and solve the related one-dimensional problem: 
we expect the characteristics displayed by the entanglement 
within the single-dot and core-shell single-dot systems to re- 
main valid (at least qualitatively) for the corresponding spheri- 
cally symmetric dots. The entanglement related to a two-wells 
one-dimensional potential corresponds instead to the one gen- 
erated within two separate spherically symmetric quantum 
dots. This system in fact respects the topology of our one- 
dimensional model. As we will see, the analysis of the re- 
lationship between the many-body wave-function characteris- 
tics and the features of the entanglement supports this choice. 

The Hamiltonian describing our system is 
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+ U{xi,X2), 



(1) 



where we have used (effective) atomic units and U{xi,X2) — 
d{xi ~ X2) models the Coulomb repulsion between two elec- 
trons in one dimension.'^ The nanostructure is modeled using 
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Dot type 


R and p range 


Potential 
(harder for larger p) 


Single doti4.i5.i6 


12 < < 2 

ZU J; it; p J; 4 


Single well 


Core-shelli^ 


9 < ^ < 16;p > 2 
16<i?<30;p>7 


Well within a well 


Double dotSiiii^ 


i? < 7; p > 1 


Double well 



TABLE I: Table showing the main R and p parameter ranges corre- 
sponding to different types of nanostructures and potentials. 



the two-center power-exponential potentialii V{x) given by 
V{x) = -K, {exp[-(|a; + d\/RY] 



exp[ 



'{\x~d\/RY]} 
(2) 

with the two centers symmetric in respect to the origin and 
situated at ±d. This allows us to study a broad class of con- 
finement potentials with different shape, size, softness and 
smoothness of the nanostructure boundaries. 

By varying its parameters, this potential can describe sin- 
gle and double quantum dot structures — to model for exam- 
ple gate-defined or self-assembled QDs — as well as core-shell 
quantum dots (such as the ones synthesized using colloidal as- 
sembling techniques'^^) or self assembled QDs with composi- 
tional modulation (see Table|I]l. 

Fig. [1] highlights how the potential varies with the value 
of the parameters, and how its flexibility allows to consider 
'soft' or 'hard' potentials as well as core-shells with different 
proportions between the two material components (compare 
for example center and bottom panel for p — 200). 

For fixed potential well depth Vb and potential range R, the 
parameter p in Eq. ^ characterizes the softness (hardness) of 
the confinement potential. If p = 1 the potential is triangular- 
like; for 2 < p < 10 the confinement evolves from a 'soft' 
Gaussian-like potential toward wells with a flat bottom and 
steep sides ('hard' potential). By increasing p even further 
(p « 200) we obtain a potential practically indistinguishable 
from rectangular-like wells. 

Each panel in Fig. [T] corresponds to a different potential 
range R, whose value determines the transition between dif- 
ferent geometries. For each R we plot the potential for three 
different values of p {p — 2,p — 7, and p = 200). 

In the case of soft confinement potential (p = 2), 
with increasing R the resulting two-center potential changes 
smoothly from two separated wells to a single potential well. 
The smoothness of this potential can simulate for example the 
controlled, incremental doping of a semiconductor structure. 
For p — 7 we obtain a steeper confinement potential, but, 
more importantly, we obtain a 'core-shell' structure for inter- 
mediate values of R. Notice that the fact that the potential is 
still relatively smooth can simulate the experimental situation 
in which there is a strong intermixing between the core and 
shell materials at their interface. For p = 200 the confine- 
ment potential is rectangular-like and the structure interme- 
diate between double and single well may correspond to col- 
loidal core-shell QD-nanostructures. In this case, we deal with 
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FIG. 1: Confinement potential profile V{x) as a function of x and 
for different values of R. Solid lines correspond to the soft confine- 
ment potential characterized by p = 2, dashed lines correspond to 
a rectangular-like confinement potential (p — 200) and dotted lines 
correspond to an intermediate confinement potential (p = 7) . Other 
parameters are Vb — 10 Hartree and d = 8 a. 



a compound QD nanostructure, which consists of the small 
inner QD embedded in a larger outer shell with no relevant 
intermixing between the different materials.'^ 

To analyze the influence of the geometry and softness 
of the confinement potential on the entanglement between 
two electrons trapped within the nanostructure, we solve the 
Schrodinger equation 



(3) 



numerically by 'exact' diagonalization. To this purpose, we 
express the wave function as a linear combination of single 
particle basis functions 

'i'i{xi,X2) = ^^aj^j^.inj^{xi;uj)nj^{x2;i^), (4) 



where {nj.{x;u!)} are the eigenfunctions of the one- 
dimensional harmonic oscillator with angular frequency uj. To 
solve Eq. (O we truncate the expansion in Eq. (|4|i by consid- 
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FIG. 2: Ground state energy E (upper panel), Coulomb repulsion 
< U > (medium panel), and linear entropy of the reduced density 
matrix L (lower panel) versus confinement potential range 7?, for 
different values of p (as labeled). The results for < (7 > and L 
corresponding to small values of R are shown in the insets of the 
middle and lower panel, respectively. 



ering only terms with 1 < < iV, / = 1, 2, where N is the 
single-particle basis size to be used in the calculation. 

We have performed our calculations for fixed Vq = 10 (ef- 
fective) Hartree and d = 8a, where a is the (effective) Bohr ra- 
dius. We are interested in the system ground state and we find 
that, using huj — 0.25 (effective) Hartree, a single-particle ba- 
sis size of = 50 is large enough to achieve convergence.^" 

For simplicity in the following we will refer to the system 
ground state and to its eigenvalue as '${xi,X2) and E respec- 
tively. 



III. GROUND STATE ENERGY AND COULOMB 
REPULSION 

To understand the entanglement properties of these sys- 
tems, it will be helpful to consider first the behavior of the 
ground state energy E and of its Coulomb repulsion compo- 
nent <U >. 



Figures 12 a) and|2b) show the results for E and < U >, 
respectively, as a function of the confinement potential range 
R and different values of p. For p = 2 ('soft' confinement 
potential) and large R the structure is a single well of depth 
2Vb exp(— d^/i?^). The two electrons are localized at the bot- 
tom of the well with an energy E « —AVoexp{—(P/R^). 
With decreasing R the energy increases as the well narrows 
and its depth diminishes. For i? w 10 the structure starts to 
split into two wells of depth of approximately Vq, centered at 
±d and separated by a shallow potential barrier. The ground 
state energy is now E « — 2Vb. The formation of the baiTier 
between the wells produces a rapid decrease of the Coulomb 
interaction down to a negligible value. The height of the bar- 
rier increases with decreasing R, while the width of each well 
decreases maintaining a depth of Vq. This reflects in a slight 
increase of E. For i? « 1 the two wells become so narrow 
that there is a rapid increase of the system kinetic energy and 
correspondingly a sharp variation in the ground state energy 
derivative: E now rapidly increases to reach its maximum 
value. As R decreases the electronic wave-function starts to 
spread out of the wells. For R < 0.1 this determines a slight 
increase of the Coulomb interaction (see inset in Fig.|3b)). 

For p = 7 and i? « 30 the confinement potential passes 
from a single well to a core-shell QD nanostructure with a thin 
outer shell and a smooth transition between core and shell. 
As R decreases, the external layer becomes thicker and more 
sharply defined, while the inner core diameter decreases in 
size. The two electrons are localized at the bottom of the inner 
well with a corresponding ground state energy of _E « — 4Vo. 
The Coulomb repulsion increases with decreasing inner core 
width. For 9 < i? < 15 the width as well as the depth of the 
inner core decreases and correspondingly both E and < U > 
rapidly increase. At i? w 10 the electronic wave function start 
to spread into the outer shell and for R < 9 the structure turns 
into a double well with a depth of Vq and an interwell barrier 
whose height increases with decreasing R. As the electrons 
localize in different wells, the Coulomb interaction suddenly 
drops to a negligible value. For i? < 1 the two separated 
wells become so naiTow that the system kinetic energy rapidly 
increases, similarly to p = 2. For i? < 0.3 the spreading of the 
electronic wave-function outside the wells determines a slight 
increase of the Coulomb interaction (see inset in Fig.|2b)). 

In the case of a rectangular-like confinement potential {p = 
200), the type of nanostructures encountered for different R 
are similar to the p = 7 case. However the increased steepness 
of the potential walls and the abrupt transition between core 
and outer shell and between the nanostructure and the sur- 
rounding material induce sharper transitions for decreasing R 
in both the ground state energy and the Coulomb repulsion. 
The energy is basically constant for R> 10 and 2 < i? < 8. 
For i? « 9 the Coulomb interaction rapidly reaches its maxi- 
mum as the inner core becomes very narrow but the electrons 
are still localized within it. Then, as the electrons delocal- 
ize into the outer shell, the Coulomb repulsion suddenly de- 
creases to a negligible value. It will slightly increase again for 
R < 0.3 as the electronic wavefunction significantly spreads 
out of the wells. 

We notice that, due to the chosen parameters, the Coulomb 
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repulsion represents at any time a very small fraction of the 
total energy. However we will show that it will play a crucial 
role for the understanding of the behavior of the entanglement. 



IV. ENTANGLEMENT AND LINEAR ENTROPY 

Bipartite entanglement is nowadays well-stated for distin- 
guishable two-component quantum systems. In particular the 
state representing an entangled system cannot be factorized 
into a product of independent states describing its parts. 

Difficulties appear in classifying and quantifying the entan- 
glement of a system composed by indistinguishable particles. 
This is due to the requirement for the antisymmetrization or 
symmetrization under particle exchange of the wavefunction 
describing indistinguishable fermions or bosons, respectively. 
In Ref. |21 22] it has been shown that the unavoidable corre- 
lations due to the particle-exchange symmetry can be related 
(for fermions) to the Slater rank of the state. This corresponds 
to the minimum number of Slater determinants in which the 
state can be expanded. For two indistinguishable fermions the 
minimum possible Slater rank is one, and the related entan- 
glement corresponds to the unavoidable antisymmetrization 
of the wave-function. This entanglement cannot be used as a 
resource for quantum-information processing. 

In the system of two indistinguishable fermions we are con- 
sidering, the entanglement is distributed over the spin and spa- 
tial degrees of freedom. In this paper we are interested in the 
entanglement of the ground state: the related many-body wave 
function can be factorized into spin and spatial components 
and the entanglement of the two parts treated separately. 

In particular the spin component of the ground state is a 
singlet state and thus always maximally entangled: if we con- 
sider then the wave function corresponding to the minimum 
entanglement content for this system, (I){xi)(f){x2){\ tii2> 
^1 iiT2>), we see that the (constant) entanglement embed- 
ded in the spin degrees of freedom corresponds indeed to that 
minimum entanglement stemming from the antisymmetry re- 
quirements. Hence, as the spin entanglement is constant, in 
the following we will focus on the entanglement generated by 
the spatial degrees of freedom only.'" 

We wish to study how the spatial entanglement changes 
as the geometry of the confinement potential is modified. 
We will quantify the spatial entanglement by using the lin- 
ear entropy L of the one-particle reduced density matrix, 
which is a useful entanglement measure for a two-fermion 
system with a very large number of degrees of freedom>i2i2^ L 
has been shown to have a behavior very similar to the Von 
Neumann entropy^'* when quantifying the particle-particle 
entanglement; '"•^■^ at the same time L is much easier to cal- 
culate, especially for a system with a very large number of 
degrees of freedom. 

The Unear entropy of the reduced density matrix is given by 



In the continuous case the reduced density matrix is given by 
with 

Pjgd(xi,a;2) = J Pred{xi,X3)pred{x3,X2)dX3 

and 

Trpia J Pl&[x,x)dx. 

The linear entropy measures the entanglement in a pure state 
by giving an indication of the number and spread of terms in 
the Schmidt decomposition of the state. 

The numerical results for the entanglement in respect to R 
are displayed in Fig. |2|c). L presents some general charac- 
teristics: a flat region with L = 0.5 for small values of R, 
followed by a sharp drop to its minimum value at interme- 
diate R, and by a partial recovery of the entanglement as R 
increases further. 

Let us consider in more detail the case of a 'soft' confine- 
ment potential (p = 2). For R > 18 the entanglement is 
lower, while the Coulomb interaction is higher, than for p = 7 
and p = 200. This is due to the shape of the p = 2 potential 
(Fig. [Tfc)) which is narrower toward the bottom of the well, 
where the electrons are confined. For the range 10 < i? < 12 
as the single QD splits into two separate potential wells the en- 
tanglement increases to its highest value and the Coulomb in- 
teraction quickly diminishes. For smaller R the overall struc- 
ture corresponds to two separated QDs and the entanglement 
does not change. The Coulomb interaction is negligible in 
this range. For R < 0.1 as the Coulomb interaction slightly 
increases, the entanglement decreases up to L « 0.3. 

For p — 7 and large R both entanglement and Coulomb 
interaction are intermediate between the cases of p = 2 and 
p — 200, L decreasing and < U > increasing with decreasing 
i?. At i? « 17.5 the curves describing < U > forp — 2 and 
p — 7 cross and so do the curves describing the entanglement: 
from Figs.|2b) and|2|c) we note that all the crossings between 
the curves describing the Coulomb repulsion correspond to 
crossings between the entanglement entropy curves, showing 
the strong correlation between the behavior of the spatial en- 
tanglement and the Coulomb interaction. For smaller R, as 
long as the inter-particle interaction increases, the entangle- 
ment decreases and drops to its minimum value, correspond- 
ing to the maximum of < U >. For 8 < i? < 9, as the 
Coulomb repulsion drops to a negligible value, the entangle- 
ment springs to its maximum value L — 0.5, the same for 
all values of pi^ L remains constant for 0.3 ^ R ^ 8, but 
for R < 0.3, each of the separated wells becomes so nar- 
row that the electronic wave function spreads out of the wells, 
so that the entanglement decreases again toward a value of 
L = 0.3. An increase of the Coulomb interaction up to 
< U >= 12 X 10^"* Hartree accompanies this last variation 
of the entanglement. 

For p = 200 the entanglement behavior is similar to the 
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previous case. However the rectangular-like confinement po- 
tential allows for a larger modulation of the electronic wave- 
function at large values of R, so in this region the entangle- 
ment is higher (and the Coulomb repulsion lower) than for 
p ~ 2 and p = 7 . For i? « 8 the inner well becomes so nar- 
row and the electrons are so confined that the entanglement 
drops to approximately zero. The minimum of the entangle- 
ment corresponds to the maximum Coulomb repulsion. As the 
Coulomb repulsion drops to zero, the entanglement increases 
sharply to its highest value, while the QD nanostructure (and 
the electronic wave-function) splits into separate wells. For 
R < 0.3 the entanglement decreases again similarly to the 
case of p = 7. Notice that the entanglement reaches a very 
similar value for Rk, {) and R « 30. 



V. THE MANY-BODY WAVE FUNCTION 

In order to explain the behavior of the entanglement for the 
different strengths of the confinement potential and values of 
the confinement potential range R, we will discuss the prop- 
erties of the system many-body wave function. 

Let us consider first p — 200 and the two extreme cases 
of maximum (L = 0.5) and minimum (L k, 0) entan- 
glement. The wave function corresponding Xo L k, Q 
(R — 8.35) is shown in Fig. [3] upper panel. Here the 
electrons are confined within the very narrow core (see the 
upper panel inset, which shows the shape of the potential). 
Due to this very strong-confinement regime, even though the 
Coulomb repulsion is maximum, the spatial part of the wave- 
function approaches the 'non-interacting' uncorrelated limit 

f{xi,X2) = 4){xi)(t>{x2), a product State which corresponds 
to no spatial entanglement. This is clearly shown by the 
comparison between the upper panel and the lower panel, 
which shows the Gaussian factorized state ^/.g(2;i: 2^2) — 
exp(— 2.T^) cxp(— 2.T2). 

At the opposite side of the entanglement spectrum, we 
find the system characterized by i? = 3.6 which corre- 
sponds to the maximum entanglement value L — Q.b. In 
this case the nanostructure is composed by two relatively nar- 
row and well-separate wells: each well would strongly con- 
fine the particles, but having more than one well provides 
an additional degree of freedom to the system. In this case 
even the weakest Coulomb interaction will then be able to 
spatially correlate the particles in a state with the structure 
^t{xi,X2) = ^[xi - d)^{x2 +d) + + d)^{x2 - d). 
This is the spatial equivalent of a spin triplet, and as such 
contains the same entanglement (L = 0.5). Fig. |4] shows 
the actual many-body wave function for R — 3.6 (upper 
panel) and, for comparison, the Gaussian 'triplet-type' state 
■^t.G{xi,X2) = exp(-(xi - 8)2/2) exp(-(x2 + 8)^/2) + 
exp(-(x2 - 8)V2) exp(-(a;i + Sf /2) (lower panel). 

The other values of the entanglement correspond to wave- 
functions intermediate between the two described. In Fig. |5] 
we present the wave-functions for p = 200 and six values 
or R. Each wave function is a three-dimensional plot in the 
coordinates (a;i, a;2, y = 4'{xi^X2))- We plot the projection of 
the wave functions on the plane (xi, y). The wave function is 




FIG. 3: Upper panel: Two-electron wavefunctions '^{x\,X2) for 
R — 8.35 andp = 200. Inset: system confinement potential. Lower 
panel: factorized function 'I'/_G(a::i, X2) = exp( — 2a;i) exp(— 2x1) 



drawn against the confining potential shape (dotted Une) and, 
for each panel, the full confinement potential is presented in 
an inset. 

The values of R have been taken as follows: R ~ 0.3 
(Fig. 121 a)) corresponds to the range where the two separated 
wells are extremely naiTow. As a consequence the electronic 
wave function spills outside and spreads over a range much 
larger than the well boundaries. For this value of R and 
within the numerical limits of our calculation, electrons are 
close to ionization (i? « Hartree). The spreading of the 
wavefunction affects the entanglement properties of the sys- 
tem: the electrons are less localized, so that their spatial entan- 
glement decreases while their Coulomb interaction increases. 
R — 3.6 (Fig. |3b)) has been discussed before: the whole 
wave-function is localized close to the bottom of the two wells 
in a 'triplet-type' state. Spatial entanglement is now maxi- 
mum as the measure of one particle in a well would imply 
certain knowledge that the other particle is in the other well. 
The description is similar for R = 7.8 (Fig.|5fc)), which cor- 
responds to a state close to the edge of the region with constant 
L = 0.5. Though the wave function shape is quite different 
from the case R = 3.6, its overall structure remains the same 
and so does its entanglement. In this respect it is interest- 
ing to notice that a wide variety of nanostructures (see e.g. 
Fig. [Ha)) or insets in Fig.|3b) and (c)) would give rise to the 
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FTG. 4: Upper panel: Two-electron wavefunctions 'I'(a;i,a;2) for 
R = 3.6 and p — 200. Inset: system confinement potential. 
Lower panel: 'triplet-type' function ^Pt.GCa;!, 2:2) — exp( — (a::i — 
8)V2)exp(-(x2 + 8)V2) + exp(-(a;2 - 8) V2) exp(-(xi + 
8)V2). 



same maximum entanglement. Surprisingly in this case the 
entanglement loses its strong sensitivity to the wave function 
details (see e.g. discussion in Ref. flU\) and the local shape 
of the potential does not affect the entanglement: Fig. |2{c) 
shows that the same value of L is obtained for 0.5 ^ i? < 8 
and 2 < p < 200. Maximum entanglement seems indeed 
to be a very robust feature and to depend only on the topo- 
logical feature of the wave function to be separable into two 
not-interconnected regions of space. 

R — 8.35 corresponds, as mentioned, to the minimum 
value L — 0.3 x lO^'^: both electrons are strongly localized in 
the (same) core well and the spatial part of the wave function 
becomes close to a factorized form. We notice that the min- 
imum of L is very sharp (see Fig. |2|c)): in contrast to max- 
imum entanglement, its minimum value can be achieved for 
very specific parameters and local confining potential shape 
only. As soon as the particles are less strongly confined, the 
system responds to Coulomb repulsion by increasing spatial 
correlation, which increases this type of entanglement. This 
explains why the minimum value of L for the softer potentials 
characterized by p = 2 and p = 7 is actually higher than for 
p = 200. 

As R increases and the core well widens, the Coulomb re- 
pulsion forces the electrons apart and correlates their posi- 
tion further: this is evident when comparing Fig. |5jd) with 
Fig. |3e) and (f). Due to the enhanced correlation the entan- 
glement increases. 

We notice that, though the confining potential is very differ- 
ent, the wave-function (and its entanglement content) is very 
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FIG. 5: Two-electron wavefunctions over the confinement potential 
wells for six values of R and p — 200, as labeled. For each system 
{R,p), the full confinement potential shape is shown in an inset. 



similar for R = 0.3 and R = 30. Again this shows that a 
similar entanglement within two electrons can be engineered 
by using very different types of nanostructures. 

The behavior of the wave functions for the softer poten- 
tials determined by p = 2 and p — 7 is similar. For p = 7 
the lesser strength of the wave-function confinement trans- 
lates into a less pronounced minimum in the entanglement; 
for p = 2 the transition between single and double well is 
even smoother, and the ground state wave-function does not 
change significantly for 15 < i? < 30 (not shown). 

We underline that both maximum and minimum entangle- 
ment correspond to a system which is strongly localized and 
strongly confined. The main difference between the two ex- 
tremes is that to achieve the maximum entanglement some 
degrees of freedom which allow for correlation are necessary 
(more than one well in this case). 



VI. ENTANGLEMENT INDICATORS 



A. Probability density at the origin 

We will now show that the entanglement behavior can be 
inferred by considering the wavefunction behavior along two 
specific directions, xi — X2 and xi — —X2- 

|^'(a;i,a::2 = a;i)|^ represents the probability density of 
finding both electrons at the same point, while |^(a;i, X2 — 
— a;i)p represents the probability density of finding the elec- 
trons at two different points which are symmetric in respect to 
the y-axis. In Fig.|6]we plot \'i>{xi,X2 — xi)\'^ (dotted line) 
and \ 'i'{xi,X2 = —xi)\'^ (solid line) forp = 200 and the same 
sequence of R values as in Fig.|5] 

Low spatial entanglement implies that there is a high prob- 
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ability of finding both particles at the same position x: this 
would in fact imply that Coulomb repulsion is unable to in- 
duce spatial correlation and the system is close to a factorized 
state. Vice-versa we expect \'i>{xi,X2 = xi)\'^ — (and, in 
particular, |5'(0, 0)p = 0) for a maximally entangled state, as 
confirmed by Fig.|6lb) and (c). 

High spatial entanglement will imply that finding one par- 
ticle at X will inform us that the other particle has a high prob- 
ability to be at a different, but specific, location, i.e., by the 
symmetry of our confinement potential, at —x. We expect 
then |^(a;i,a;2 ~ —xi)\'^ to present peaks symmetric in re- 
spect to xi = when the spatial entanglement is non-zero, as 
shown in Fig.|6] 

For R corresponding to minimum entanglement, 
|*(xi,a;2 = -a;i)|2 « |*(a;i,a;2 = a;i)|2 and |*(0,0)|2 
reaches its maximum (Fig. |6td)). The lack of differentiation 
between these two 'orthogonal' type of correlations (having 
the same probability of finding a particle at the same position 
or in a position symmetric in respect to the origin) indicates 
in fact that little information on the other particle can be 
gained by measuring the position of one of the particles. By 
looking at Fig. |6jd) we notice once more that, for minimum 
entanglement the wave function assumes a form close to the 
factorized form ^(xi, a;2 = xi) « $(xi)<i>(a;2). 

In general the value of the probability density \'^{xi, X2)\'^ 
at the highly symmetric point xi = X2 = (common to both 
the directions considered) will increase for decreasing entan- 
glement: let us consider the values of L and |5'(0, 0)^ for 
R = 15 and the different values of p. From Fig. |2c) we see 
that L{p = 7, i? = 15) < L{p = 2, i? = 15) < L{p = 
200, R = 15). By comparing Fig.lTJa), (c) and Fig.|6je) we 
find that instead |^'(0,0;p = 7,R = 15)\^ > |*(0,0;p = 
2,R = 15)|2 > |^'(0,0;p = 200, i? = 15)1^. If we now 
consider R — 30, the p = 2 and p = 7 curves for the entan- 
glement have crossed, and indeed we find that L{p — 2,R ^ 
30) < L{p = 7,R = 30) < L{p = 200, R = 30), and, 
from Fig. Hb), (d) and Fig.lSf), |*(0, 0;p = 2, i? = 30)1^ > 
|*(0,0;p = 7,R^ 30)|2 > |*(0,0;p = 200, i? = 30)1^. 
These findings confirm that the value of the probability den- 
sity at a single but highly symmetric point contains very rel- 
evant information on the value of the overall system entan- 
glement. False indication of entanglement might occur when 
considering mixed states, but similar problems would be en- 
countered when using pure-state entanglement measures such 
as L or the Von Neumann entropy, and even criteria designed 
for mixed states will not always detect bipartite entanglement 
between systems with large numbers of degrees of freedom.— 



B. Coulomb interaction and effects of long-range repulsion 

So far we have modeled the Coulomb repulsion as a con- 
tact interaction. This implies that the Coulomb repulsion is 
non zero only if the wave function allows for both particles 
to be at the same place. The higher the probability of having 
particles at the same position, the higher the Coulomb repul- 
sion and the lower the entanglement. This counterintuitive 
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FIG. 6: \'i'{xi,X2 = xif) (dotted line) and |>I>(a;i,a;2 = 
(solid line) in respect to xi for p = 200 and the same sequence of 
values of R as in Fig.|5] For clarity a zoom of the wave-function cuts 
around a::i = are plotted in the inset of panel (d). 
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'inverse' relationship between Coulomb repulsion and entan- 
glement behaves as a very good entanglement indicator, as 
can be observed by comparing Fig. |2b) and (c), where the 
maxima (minima) of the Coulomb interaction corresponds to 
the minima (maxima) of the entanglement. We may wonder 
though if this is an artefact of the contact-type of interaction 
used. 

In Fig. [8] we present the calculations done using the long 
range Coulomb repulsion 



Ulr{xi,X2) 
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y^l + ixi -X2f 
and a basis size N = 30, which allows for convergency for 
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FIG. 8: Coulomb repulsion < 17 > and entanglement entropy L 
versus potential range R for the long-range Coulomb repulsion[5]and 
the Gaussian-type confinement potential p = 2 

the range of values of R shown. We choose p = 2 which, de- 
scribing the softest potential among the set of p- values chosen, 
would allow for the biggest modification of the corresponding 
many-body wave-functions, and hence of the entanglement. 

Fig. [8] shows first of all that the main characteristics of 
the entanglement found with the contact interaction are con- 
firmed: the entanglement entropy presents a plateau with 
L = 0.5 for i? < 10, a rapid decrease for intermediate values 
of R reaching a minimum for i? « 17, and a slow increase of 
the entanglement as R is increased further. The higher value 
of the minimum, in respect to the results obtained using the 
contact-type of interaction, witnesses the increase of spatial 
coiTelations due to the long range nature of Eq. (|5]l. 

Most importantly though, the results obtained using a long 
range interaction still show the same 'inverse' correlation be- 
tween Coulomb repulsion and entanglement: this indicates 
that Coulomb repulsion between particles, is a good indica- 
tor for spatial entanglement. 



VII. POTENTIAL QUANTUM PHASE TRANSITION 

A point of nonanalyticity in the ground state energy 
of a quantum system is associated with a quantum phase 
transition.-^ In Ref. |28] it was shown that in such a transi- 
tion a nonanalyticity of the entanglement measure would be 
associated to the nonanalyticity of the ground state energy. 

In the system we are considering, for increasing p, as the 
confining potential becomes harder and harder, a discontinuity 
in the derivative in respect to the potential range of the ground 
state energy and of the entanglement measure, dE/dR and 
dL/dR, seems to appear. This discontinuity underlines the 
transition between minimum and maximum values of the en- 
tanglement (see i? « 8 forp = 200). 

A similar pattern for the entanglement was observed in the 
quantum phase transition for two electrons close to the ion- 
ization point of a single QD.^ In that case the transition was 
between bound and unbound (resonance) states, while in our 
case a transition seems to occur between two different sets of 
bound states. 

In the system we are considering, the transition is triggered 



by a shape-change in the confining potential, from the poten- 
tial in the inset of Fig. |3d) to the one in the inset of Fig. |3c). 
Due to this change, the system shifts between two very dif- 
ferent sets of ground states: the first set describes the par- 
ticles being highly confined in the naiTow core region of a 
core-shell type structure, and it is formed by wave functions 
similar to the factorized 'i/ f{xi, X2)', the second set describes 
particles confined within two separate wells and it is formed 
by the topologically different 'triplet-type' states 'i/t{xi,X2)- 
The fact that the system ground state on the left and right of 
the transition has well defined but very different properties is 
consistent with a quantum phase transition picture (see e.g. 
Ref. Liai). 

We note that as R decreases and the transition is ap- 
proached, the energy difference between the system ground 
state, as bounded within the inner well, and the ground state 
that the system would have if the inner well would be absent 
(and which would coiTespond to a double well with a baiTier 
of vanishing width^") decreases as well. This energy distance 
between the two relevant sets of bound states reaches a mini- 
mum at the transition point. 

The transition between these two confinement potential 
shapes could be induced experimentally by changing the po- 
larity of a gate positioned over the core-to-barrier region 
within a gate-defined QD. Recent studies have shown the pos- 
sibility of engineering gate-defined QD confinement poten- 
tials with shapes ranging from Gaussian, to rectangular-type 
potentials?^ so this type of device should allow the transi- 
tion between minimum and maximum entanglement to be ex- 
plored even in (or at least close to) the large p hmit. 



VIII. CONCLUSIONS 

We have studied the entanglement of two interacting elec- 
trons confined within single, core-shell and double quantum 
dots. The confinement potential has been parametrized by a 
two-center power-exponential potential, which has allowed us 
to investigate quantum dots described by either hard or soft 
potentials, with different ranges, including the effects of the 
transition between the different types of structures. The calcu- 
lation has been done by direct diagonalization of the Hamil- 
tonian including the Coulomb interaction between the elec- 
trons. By varying the confinement potential as a function of 
dot shape, range of confinement potential — which determines 
the QD size — and the strength of the confining potential we 
showed that it is possible to induce fast and large variations 
of the entanglement between the two electrons. This prop- 
erty may be used to design nanostructures — and nanostruc- 
ture modulations via external fields — according to the level of 
entanglement required by a specific application. 

We have studied in detail the relationship between Coulomb 
repulsion and spatial entanglement and shown that they dis- 
play a counterintuitive 'inverse' correlation: due to Coulomb 
repulsion, electrons tend to correlate their position to mini- 
mize their interaction, which implies minimizing the probabil- 
ity that electrons could be found at the same position. From 
the entanglement point of view this correlation means enhanc- 
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ing the probability that if one electron is measured at a cer- 
tain position the other will be in a different but correlated 
region, thus enhancing the spatial entanglement between the 
two particles. We note however that, if the Coulomb inter- 
action is switched off, the spatial entanglement vanishes (see 
also Ref. ll32ll '). This is in contrast to the behavior of the 'lo- 
cal' (or 'site') entanglement which characterizes the Hubbard 
model, where zero Coulomb interaction corresponds to maxi- 
mum entanglement. 

We have analyzed the many -body wave function and in par- 
ticular the correlations between the particle probability den- 
sity along some specific directions and the entanglement. We 
have then proposed the value of the particle probability den- 
sity at a single (but highly symmetric point) as an indicator of 
entanglement. 

We have identified a potential quantum phase transition be- 
tween minimally and maximally entangled states within our 
system. This transition is triggered by a change in the poten- 



tial shape which induces a topological change in the many- 
body wave function, from an almost factorized to a 'triplet- 
type' wave-function form. Further investigation of this in- 
triguing phenomena will be pursued in future work. 

Systems of electrons confined in quantum dots have 
been proposed as tools for performing quantum informa- 
tion/computation tasks; it is then of great importance to un- 
derstand how the entanglement between these particles can be 
engineered and tailored. Our work provide a systematic study 
in this direction. 
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